Single-Particle Momentum Distribution of an Efimov trimer 
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Experimental progress in the study of strongly interacting ultracold atoms has recently allowed the 
observation of Efimov trimers. We study theoretically a non-conventional observable for these trimer 
states, that may be accessed experimentally, the momentum distribution n(k) of the constitutive 
bosonic particles. The large momentum part of the distribution is particularly intriguing: In addition 
to the expected tail associated to contact interactions, it exhibits a subleading tail which 
is a hall-mark of Efimov physics and leads to a breakdown of a previously proposed expression of 
the energy as a functional of the momentum distribution. 

PACS numbers: 31.15.ac,03.75.Hh 



C/2 



o 



(N 

> 

in 

rn 
O 



X 



I. INTRODUCTION 



Experiments with ultracold atoms have now entered 
the regime of strong interactions, thanks to the possibil- 
ity to manipulate the s-wave scattering length a between 
cold atoms with a magnetically induced Feshbach reso- 
nance [H, Q] ■ This has led to a revolution in the study 
of the few-body problem, as one can now have in a con- 
trollable way a scattering length much larger (in abso- 
lute value) than the range h (and the effective range) 
of the interaction potential. In particular, this has al- 
lowed to confirm experimentallyll, 0] the existence of 
the long-searched Efimov effect p-uf: As shown by Efi- 
mov in the early 1970's, three particles interacting via a 
short range potential with an infinite scattering length 
may exhibit an infinite number of trimer states with a 
geometric spectrum. The existence of an infinite num- 
ber of bound states is usual, even at the two-body level, 
for long range interactions, but it is quite intriguing for 
short range interaction potentials. This Efimov effect 
takes place for three (same spin state) bosons Q, but it 
is more general, it also occurs for example for two (same 
spin state) fermions and a third distinguishable particle 
at least 13.607 times lighter 

On the experimental side, an increasing number of ob- 
servable quantities are now at hand. For Efimov physics, 
the usual evidence of the emergence of an Efimov trimer 
state is a peak in the three-body loss rate as a func- 
tion of the scattering length Q. Now radio-frequency 
spectroscopic techniques can give a direct access to the 
trimer spectrum ^] . For strongly interacting Fermi gases 
(without Efimov effect) a very precise measurement of the 
atomic momentum distribution n(k) was performed re- 
cently, so precise that it allowed to see the large momen- 
tum tail n(k) ^ C/k^, k large but still smaller than 1/6, 
and to quantitatively extract the coefficient C whose val- 
ues were satisfactorily compared to theory ^] . The same 
conclusion holds for the few-body numerical experiment 
of |10|]. Similarly, the first order coherence function g^^^ 



of the atomic field over a distance r, a quantity measured 
for bosonic cold atoms [ll| but not yet for fermionic cold 
atoms, is related to the Fourier transform of ri(k) and is 
sensitive to the 1/k^ tail by a contribution that is non- 
differentiable with respect to the vector r in r = p^ . 
and that ap peared in the many-body numerical experi- 
ment of [l3| . 

The occurrence of the l//c^ tail in n(k) is a direct con- 
sequence of two-body physics, that is of the binary zero- 
range interaction between two particles, and it holds in 
all spatial dimensions: According to Schrodinger's equa- 
tion for the zero-energy scattering state 0(r) of two par- 
ticles of relative coordinates r, Ar(/)(r) oc 5{v) for a con- 
tact (regularized Dirac delta) interaction [14|, so that in 
Fourier space ^(k) oc l/Zc^ and n(k) cx |^(k)p scales as 
1/k'^. On the contrary, the coefficient C, called contact, 
depends on the many-body properties, and can be related 
to the derivative of the gas mean energy (or mean free 
energy at non-zero temperature) with respect to the scat- 
tering length, as was shown first for bosons in one dimen- 
sion Il5|, then for spin 1/2 fermions in three dimensions 
[T6l - ll8l |. and for bosonic or fermionic, three dimensional 
or bidimensional systems, in 19]. 

In this paper, we anticipate that experimentally, it may 
be possible to measure with high precision the atomic mo- 
mentum distribution in systems subjected to the Efimov 
effect, for exarnple in a Bose gas with a large scatter- 
ing length |2C)| - [23| . To be specific, we consider in the 
center of mass frame the Efimov trimer states for three 
bosons interacting with infinite scattering length. After 
recalling the expression of the three-body wavefunction 
in section ini we obtain the expression of the momentum 
distribution in terms of integrals over a single momentum 
vector in section |TTT1 see Eqs. (|18ll9l20l21l22l) . As iUus- 
trated in section lTVl this allows to perform a very precise 
numerical evaluation of the momentum distribution for 
all values of the single-particle wavevector k, and to an- 
alytically obtain the large momentum behavior of n(k): 
In addition to the expected C /k'^ term at large A:, we find 



2 



an unexpected subleading term, that is a direct and 
generic signature of Efimov physics [l^, see An- 
other, more formal, consequence of this 1/fc^ subleading 
term is that the general expression giving the energy as 
a functional of the momentum distribution rt(k), derived 
for the non-Efimovian case in [25] and extended (with 
the same form) to the Efimovian case in [26l |. turns out 
to be invalid in the Efimovian case [131 ■ We conclude in 
section IVl 



ip is related to the Jacobi coordinates r = r2 — ri and 
p = [2r3 — (ri + r2)]/-\/3. In our expression of VE*, Vcm is 
the Gaussian wavefunction of the center of mass ground 
state in the harmonic trap, normalized to unity, and tp is 
a Faddeev component of the free space trimer wavefunc- 
tion. The explicit expression of ^ is known Q: 



w ^ _ A/"v, K,„{Kay/r^ + p^) sin[so(f - a)] 



sin(2Q:) 



II. NORMALIZED WAVEFUNCTION OF AN 
EFIMOV TRIMER 

A. Three-body state in position space 

In this subsection we recall the wavefunction of an Efi- 
mov trimer and give the expression of its normalization 
constant. We consider an Efimov trimer state for three 
same-spin-state bosons of mass m interacting via a zero 
range potential with infinite scattering length. In or- 
der to avoid formal normalisability problems, it is con- 
venient to imagine that the Efimov trimer is trapped in 
an arbitrarily weak harmonic potential, that is with a 
ground state harmonic oscillator length a^o arbitrarily 
larger than the trimer size [28]. Since the center of mass 
of the system is separable in a harmonic potential, this 
fixes the normalisability problem without affecting the in- 
ternal wavefunction of the trimer in the limit Oho +oo. 
In this limit, the energy of the trimer is equal to the free 
space energy 



£^trim = KO > 0. 



(1) 



According to Efimov's asymptotic, zero-range theory 



V2 
Rt 



g 7r9/|so|gArgr(l-hso)/|so| 



(2) 



where R * > is a length known as the three-body pa- 
rameter [3l|, the quantum number q is any integer in Z, 
and the purely imaginary number sq = j|so| is such that 



where Ks„ is a Bessel function and a — atan(r/p). The 
normalization factor Af^p ensuring that lll'lp = 1 was 
not calculated in Q. To obtain its explicit expression, 
one first performs the change of variables (ri,r2,r3) — > 
(C,p, r), whose Jacobian is £'(ri, r2, r3)/D(C, p, r) = 
(— \/3/2)'^. To integrate over r and p one then intro- 
duces hyperspherical coordinates in which the wavefunc- 
tion separates; one then faces known integrals on the 
hyperradius [32] and on the hyperanglcs 33] . This leads 
to HI: 



2nl cosh(^) 



cosh( ) H smh( ; 
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cosh( — - — ^ 



3V^ 



(6) 



B. Three-body state in momentum space 



To obtain the momentum distribution for the Efimov 
trimer, we need to evaluate the Fourier transform of the 
trimer wavefunction ^ given by ^ . Rather than directly 
using ([U, we take advantage of the fact that, for contact 
interactions, the Faddeev component ip obeys the non- 
interacting Schrodinger's equation with a source term. 
With the change to Jacobi coordinates, the Laplace op- 
erator in the coordinate space of dimension nine reads 



2[Ar 



so that 



[4 



Aj^(r,p) - J(r)B(p). 



(7) 



jsol cosh(^-) = -j= smh(-^), (3) 

so that jsol = 1.00623782510... The corresponding 
three-body wavefunction may be written for Koiho — >■ 
-l-oo as 



*(ri,r2,r3) ?Acm(C) 
|2ri-(r2+r3)| 



■ |2r3-(ri+r2)| , 
i^{r.., ) 

2r2 - (rs +ri)| 



V3 



V3 



(4) 



The source term in the right hand side originates from 

the fact that 



'4}{r, p) 



B{p) 



(8) 



for a fixed p, this 1 /r divergence coming from the replace- 
ment of the interaction potential by the Bethe-Peierls 
contact condition. Taking the Fourier transform of ([7]) 
over r and p leads to 



B(K) 



(9) 



where C = (ri -I- r2 -I- r3)/3 is the center of mass po- 
sition of the three particles and the parameterization of 



where the Fourier transform is defined as B{K) = 
J (P pe~^^'''B{p). B{p) is readily obtained from ([5]) by 
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taking the limit r 0: 

Bip) = -AA,(4.)V^.sinh(^)^f^. (10) 

z p 

The Fourier transform of this expression is known, see 
relation 6.671(5) in [s^], so that 



B{K) = 



X(if2 + K2)l/2 



Ko 



(11) 



where c.c. stands for the complex conjugate. Note that 
the expression between the curly brackets simply reduces 
to 2i sin(|so|Q!) if one sets K = Kosinha. What we shall 
need is the large K behavior of B{K). Expanding (fTTj) 
in powers of n^/ K gives 



B{K)^N^'^-^[{2K/n^) 



+ 0{l/K^). (12) 



The last step to obtain the trimer state vector in mo- 
mentum space is to take the Fourier transform of (|H), 
using the appropriate Jacobi coordinates for each Fad- 
deev component, or simply by Fourier transforming the 
first Faddeev component using the coordinates (C, r, p) 
given above and by performing circular permutations on 
the particle labels. This gives 



*(ki,k2,k3)- 



kil |2k3 



?/'CM(ki+k2 

(ki+k2)| 



2^3 



ka) 

(231) + (312) 



(13) 



where the notation (ijfc) means that the indices 1,2,3 
have been replaced by k respectively. 



III. INTEGRAL EXPRESSION OF THE 
MOMENTUM DISTRIBUTION 

To obtain the momentum distribution for an Efimov 
trimer state, it remains to integrate over ks and k2 
the modulus square of the Fourier transform (|13p of the 
trimer wavefunction. In the limit KoOho ~> +oo where 
one suppresses the harmonic trapping, one can set 

|V;cM(ki + k2 + k3)P - {2^f5{\ii + k2 + ks) (14) 

so that the trimer is at rest in all what follows. Integra- 
tion over ks is then straightforward: 



r Jkz-kil V3|ki+k2| 



. ^(|k2 + ikil, ^fci) + V^dki + ik2|, ^fc2)f . (15) 



The factor 3 in the right hand side results from the fact 
that, as e.g. in [l6j . we normalize the momentum distri- 
bution n(k) to the total number of particles (rather than 
to unity): 



(27r)3 



n(k) 



(16) 



Also note that the sum of the squares of the arguments 
of is constant and equal to k\ + k"^ k2 for each 

term in the right hand side of (IT5|) . When using 
one can thus put the denominator in Q as a common 
denominator, to obtain 



n(ki) 



(Pk. 



^(#|ki 



k2|) 



S(^fci) 



B{^k2) 



(27r)3 



(43/34)(fc2+A;|+ki •k2 + K2)2 



(17) 



For simplicity, we have assumed that the normalization 
factor A/"^ is purely imaginary, so that B{K) is a real 
quantity. 

In the above writing of ?i(ki), the only "nasty" contri- 
bution is _B(\/3|ki + k2|/2); the other contributions are 
"nice" since they only depend on the moduli ki and /c2- 
Expanding the square in the numerator of (jl7p . one gets 
six terms, three squared terms and three crossed terms. 
The change of variable k2 = — (k2 + ki) allows, in one 
of the squared terms and in one of the crossed terms, to 
transform a nasty term into a nice term. What remains 
is a nasty crossed term that cannot be turned into a nice 
one; in that term, as a compromise, one performs the 
change of variable k2 — — (k2 + ki/2). We finally obtain 
the momentum distribution as the sum of four contribu- 
tions, 

n(ki) = ni(ki) -|- nii(ki) + niii(ki) + niv(ki), (18) 
with 



ni(ki) = 
nii(ki) = 
niii(ki) = 
niv(ki) = 



BH^k,) 



437 (2^)3 (fc2 + kl + ki • k2 + ^2)2 
34 fd^k2 2B\^k2) 



437 (2^)3 (fc2 + kl + ki • k2 + nlY 
34 /•d3fc2 AB{^ki)B{^k2) 



437 (2^)3 (fc2 + A:2 + ki • k2 + nlY 



(19) 
(20) 
(21) 



34 [d^k2 2B(^|2k2+ki|)B(^|2k2-ki| 



43 J (27r)3 



(kq + k2 + 1^1)^ 



(22) 



An interesting question is to know if one can go be- 
yond the integral expressions Eqs. (|19l20l21l22l) . that is if 
one can obtain an explicit expression for the momentum 
distribution, at most in terms of special functions. The 
contribution ni(ki) is straightforward to calculate: 



, V3 _BV3fci/2^ 
4^ Uy (fc? + 4^2/3)1/2- 



(23) 
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The contribution nii(ki) is also explicitly calculable by 
performing the change of variable k2 = (2/-\/3)ko sinha 
and using the identity that can be derived from contour 
integration: 



+00 



da 



e*'*" 27r sin[s(i7r — ao)] 

cosh a — coshao sinhag sinh(s7r) 



(24) 



where s is any real number and uq is any complex num- 
ber with < Im ofQ < vr. This also allows to obtain an 
explicit expression of nin(ki) if one further applies inte- 
gration by part, integrating the factor sin(|so|a). We do 
not give however the resulting expressions since, contrar- 
ily to these first three contributions to ri(ki), the contri- 
bution niv(ki) in (22) blocked our attempt to calculate 
n(ki) explicitly. For ki = however it becomes equal 
to the contribution nn. n(ki = 0) can thus be evaluated 
explicitly in terms of kq and sq: see Appendix F in p^ . 
In numerical form this gives 



n(k = 0) 



55.43379775608.. 



(25) 



IV. APPLICATIONS 

A. Numerical evaluation of n(k) at all k 

The integral expression of 7i(k) derived in section IIIII 
allows a straightforward and very precise numerical cal- 
culation of the single-particle momentum distribution 
for an infinite scattering length Efimov trimer, once all 
the doable angular integrations have been performed in 
spherical coordinates of polar axis ki. The result is 
shown for low values of k in FigH^, and for high values 
of k in FiglTJ). In particular, Fig[T}3 was constructed to 
show how n(k) approaches the asymptotic behavior ((26|) 
derived in the next subsection, that is to reveal the ex- 
istence of a sub- leading oscillating term. Note that 
the accuracy of the numerics may be tested from (j25p 
and from the explicit analytical expressions of ni (given 
in (j23p ). of Till and of nm (not given). 



B. Large momentum behavior of 7i(k) 

Starting from the integral representation 
Eas. (|19l20l21l22p . we show in the Appendix \M that 
the single-particle momentum distribution has the 
asymptotic expansion at large wavevectors: 



2|so| ln(-\/3fc/Ko) + (p 



(26) 



where we recall that the trimer energy is i^trim — 
—h^K^/m, and the quantities C, D and (p derived in the 



Appendix |X] are given by 



C/ko = 87r^sinh(|sok/2)tanh(|so|7r)/ cosh(^-^) 



+ __smh(^-) - 

53.09722846003081 

-89.26260 

-0.8727976. 



47r , Jsok> 



(27) 

(28) 
(29) 
(30) 



The crucial point is that I? 7^ 0: Due to the Efimov 
effect, the momentum distribution has a slowly decaying 
0(l/fc^) oscillatory subleading tail. 



Breakdown of the usual energy-momentum 
distribution relation 



In [26| it was proposed that the expression of the en- 
ergy as a functional of the momentum distribution, de- 
rived in [2^ for equal mass spin 1/2 fermions, also holds 
for bosons (apart from the appropriate change of numer- 
ical factors). In the present case of a free space Efimov 
trimer at rest with an infinite scattering length, the en- 
ergy formula of ^] reduces to 



(fk ti'e 



(27r)3 2m 



n(k) 



C_ 



(31) 



We have however put a question mark, because the 
asymptotic expansion (j26p implies that the integral in 
(Plj) is not well-defined: After the change of variables 
X — ln(\/3fc/Ko), the integrand behaves for x — -l-oo as a 
linear superposition of e^*'*"'^ and g-^A'^ol^^ ^jj^g^^ jg 3,3 a 
periodic function of x oscillating around zero. This was 
overlooked in [2^ . 

At first sight, however, this does not look too serious: 
One often argues, when one faces the integral of such an 
oscillating function of zero mean, that the oscillations at 
infinity simply average to zero. More precisely, let us 
define the cut-off dependent energy functional 



k<A 



(27r)3 2m 



n(k) 



C_ 

fc4 



(32) 



where the integration is limited to wavevectors k of mod- 
ulus less than the cut-off. For A — > +00, E{K) is asymp- 
totically an oscillating function of the logarithm of A, 
oscillating around a mean value E. The naive expec- 
tation would be that the trimer energy -Etrim equals E. 
This naive expectation is equivalent to the usual trick 
used to regularize oscillating integrals, consisting here in 
introducing a convergence factor g^-'n^r^(^Jlk/^^Q) |-j^g 
tegral without momentum cut-off and then taking the 
limit 7^ — > 0^: 



lim 



(27r)3 2m 



"(k) - 



E. 

(33) 
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FIG. 1: For a free space Efimov trimer at rest composed of three bosonic particles of mass m interacting via a zero-range, 
infinite scattering length potential, single-particle momentum distribution n(k) as a function of k. (a) Numerical calculation 
from the expression for n(k) appearing in Eq. (|18p . (b) Numerical calculation (solid line) and asymptotic behavior (|26|l (dashed 
line), with the horizontal axis in log scale. The unit of momentum is kq, such that the trimer energy is —f?KQ/m. 



To test the naive regularization procedure, we per- 
formed a numerical calculation of £'(A), using the result 
([TS]) to perform a very accurate numerical calculation of 
n(k). The result is shown as a solid line in FigH] We also 
developed a more direct technique allowing a numerical 
calculation of ^^(A) without the knowledge of n(k), see 
Appendix |B] The corresponding results are represented 
as -|- symbols in Figl5] and are in perfect agreement with 
the solid line. As expected, E{K) is asymptotically an os- 
cillating function of the logarithm of A, oscillating around 
a mean value E. 

To formalize, we introduce an arbitrary, non-zero value 
fcniin of the momentum. We define 5n{k) = ri(k) — C/k'^ 
for k < fcmin, and for k > fcmin: 



Snik) = n(k) - 7T + 



■ cos 



2|so| \n{V3k/Ko) + (p 



C_ D_ 

(34) 

The introduction of fcmin ensures that the integral of 
k'^Sn{k) over all k converges around k = 0. The sub- 
traction of the asymptotic behavior of n(k) up to order 
0{l/k^) for k > kmin ensures that the integral of k'^Sn{k) 
over all k converges at infinity. As a consequence we get 
in the large cut-off limit 



E{A) = E- 



- sm 



[2|so| ln(V3A/Ko) + <^] + 0(l/A), 



with 



E 



(35) 



■ sm 



[2|so| ln(V3fc,nin/Ko) + ip] 







+00 ^2l4 

dk—^5n{k). (36) 



From this last equation (j36p and the numerical calcu- 
lations of n(k) first up to A; = IOOOkq and then up to 
k ~ 5500koj we get two slightly different values of E, 
which gives an estimate with an error bar [35| : 



E ~ 0.89397(3)i;tri 



(37) 



The conckision is that E (significantly) differs from -Btrim: 
The naive regularization of the energy formula proposed 
in 26] does not give the correct value of the trimer energy. 



An analytical representation of E in terms of single 
integrals can be obtained, see Appendix This analyti- 
cal calculation gives a physical explanation of the failure 
of the naive regularization: It is inconsistent to add by 
hand the regularization factor g^-^^^'^i^^/'^o) a,t the last 
stage, that is in the integrand of (j31[) . To be consistent, 
the momentum cut-off function has to be introduced at 
the level of the three-body problem. Then the subleading 
term in the momentum distribution acquires a small 
non-oscillating component, of order rj/k^, that gives a 
non-zero contribution to the integral (|5T|) for — > 0+, 



4:Vp-rilniV3k/Ko) 



does not tend to zero in 



since J^, dkk^-^e 
this limit. The resulting integral representation of E con- 
firms the numerical result and allows to evaluate E with 
a better precision: 
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1 



4cosh(|so|7r/2) 



cosh(|so|7r/2) 



|sok 



■ sinh(|so|7r/2) 



^cosh(|sok/6) 



^ ^+°°^^ f 21n(g)cos(|so|lng) ^ 16gln[(l + g2)/4] 



1 + g2 + ^4 



h41n 



1 



= 0.8939667780883 



cos [ I So I In 



1 



'^)/4] . /, I, l + q 

1 TTT-r^ TTT Sill Sn M 1 ; 

1 , ln[2(l + 92)/(g2 + 3)] 



2(92 + 3) 



1 -g2 



(38) 
(39) 




FIG. 2: Cut-ofT dependent energy £;(A) as defined in ((Ml) 
for a fi'ee space infinite scattering length Efimov trimer with 
a zero-range interaction, as a function of the logarithm of 
the momentum cut-off A. Solid line: numerical result ob- 
tained via a calculation of the momentum distribution n(k). 
Symbols -|-: direct numerical calculation of E{K) as exposed 
in the Appendix [B] Dashed sinusoidal line: asymptotic os- 
cillatory behavior of -E(A) for large A, obtained by omitting 
0(1/A) in H35p . Dashed horizontal line: mean value E around 
which E{K) oscillates at large A. The values of E obtained 
analytically (|38} and numerically (|37|) are indistinguishable 
at the scale of the figure, and clearly deviate from the dot- 
ted line giving the true energy i?trim of the trimer, exempli- 
fying the failure of a at a first sight convincing application 
of an energy formula for bosons in three dimensions. The 
unit of momentum kq is such that the true trimer energy is 

-Etrim = -ft^Ko/m. 



V. CONCLUSION 

We have calculated the single-particle momentum dis- 
tribution n(k) for the free space Efimov trimer states 
of same spin state bosons interacting via a zero-range 
potential with an infinite scattering length. The asymp- 
totic behavior of n(k) at large wavevectors, that we de- 
termined with good precision, is of particular interest: 
In addition to the C/A:^ tail expected from two-body 
physics, it has a subleading oscillating contribution, 
which is a signature of Efimov physics that one may try 
to observe experimentally. 

We obtained the analytical expression for the coeffi- 
cient C, see (|?7)) . This coefficient can also be obtained 



by a direct calculation in position space, usiiig the fact 
that C is proportional to / (fp\B{p)\^ [H, This 
result allows to calculate the trimer energy for a finite 
scattering length a, to first order in 1/a, thanks to the 
relation 11911 



dEtii-m 



(40) 



where the derivative is taken for a fixed value of the three- 
body parameter Rt- In other words, we obtained analyt- 
ically the derivative at —tt/2 of Efimov's universal func- 
tion A(^) In numerical form, it gives A'(— 7r/2) — 
2.125850069373..., which refines the previously known 
numerical estimate ~ 2.12 5]. Furthermore, the exis- 
tence of the 1/fc^ subleading term leads to a failure of 
the relation proposed in 26!] expressing the energy as a 
functional of the single-particle momentum distribution, 
which was not obvious a priori. 

We have considered here the particular case of Efimov 
trimers. The coefficient of the subleading term was 
however obtained in Appendix |X] by taking the zero- 
energy limit kq 0. We thus expect that the phe- 
nomenology of the subleading term persists, not 
only for any other three-body system subjected to Efimov 
physics (such as two fermions of mass M and a lighter 
particle of mass m with a mass ratio M/m > 13.607 
0, 1^), but also for a macroscopic Bose gas, at least for 
strong enough interactions [U [2^ so as to make the 
subleading 1/fc^ term sizeable and maybe accessible to 
measurements. More precisely, we conjectured in [l^ 
that there is a subleading oscillating term in the 
tail of the momentum distribution [as in Eq. (j26p ] for any 
A^-body quantum state subjected to the Efimov effect, 
that is obeying three-body contact conditions involving 
the three-body parameter Rt , with a coefficient D related 
to dE / d{hvRt) through a simple proportionality factor. 
Moreover, this derivative was related in J^] to a three- 
body analog of the contact, defined from the prefactor 
appearing in the three-body contact condition. 
Note: After submission of this work, (i) the expression 
of the coefficient of the 1/fc^ subleading tail of n(k) in 
terms of the derivative of the energy with respect to the 
logarithm of the three-body parameter Rt, and (ii) the 
appropriate energy formula taking into account the Efi- 
movian subleading 1/fc^ term, appeared in [36] for a Bose 
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gas with an arbitrary number of particles. 
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Appendix A: Leading and next-to-leading terms for 
n(k) at large momentum 



Here we derive the asymptotic expansion (|26|) . We 
shall take the large ki limit, or equivalently formally the 
Ko — > limit for a fixed fci. From the asymptotic be- 
havior ([T^ we see that B^(ki) involves a sum of "os- 



cillating" terms involving k^^° or , and of "non- 

oscillating" terms. We shall calculate first the resulting 
non-oscillating contribution, then the resulting oscillat- 
ing one, up to order included. 



a. Non-oscillating contribution up to 0{l/kf) 

We consider the small kq limit successively for each of 
the four components of n{ki) in ()18|) . 
Contribution I: Taking directly — > in the integral 
defining nj, replacing B{ki) by its asymptotic behavior 
(IT^ and averaging out the oscillating terms fcf gives 



the leading behavior 



(Al) 



Contribution II: In the integrand of (PU)) . we use the 

splitting 



[k 



ki • k2 + kI) 2 = fci 



^2 



ka + nlY 



(A2) 



The first term in the right hand side gives a contribution 
exactly scaling as 1/fcf. In the contribution of the second 
term in the right hand side, one may take the limit — >■ 
and replace B^{V3k2/2) by its asymptotic expression 
to get the subleading 1/fcf contribution. Performing the 
change of variable ka = fciq in the integral and averaging 
out the oscillating terms kf^^° gives 

{niiik,)) = ^ - ^|A/-^p^ + o{l/k!), (A3) 



with 



C = 3(V3/2)«| ^2B'iV3k 



:/2). (A4) 



B{kq sinh a) 
This leads to 



Kq sinh a cosh a 



(A5) 



C= 127r3(V3/2) 



3IW 

Kq 



+00 2 

da — 



(e^^o" -I- c.c.) 



cosh a 



(A6) 

where we used the fact that J\f^ — — \J\f^\'^. The resulting 
integral over a may be extended over the whole real axis 
because the integrand is an even function of a; it may 
then be evaluated by using the general result (that we 
obtained with contour integration) 



K{0,.s) 



+00 



da ■ 



2-K sinh(s6') 



cosh a + cos 6* sin6' sinh(s7r) 

(A7) 

where s is a real number and d e]0,7r[. One simply has 
to take 6' = 7r/2, s^O and s = 2\so\ respectively. We 
get 



24^/^^^ sinh^(|go|7r/2) „ 
no [2 J cosh(|.ok) l-^'^l • 



(AS) 



This, together with ([6]), leads to the explicit expression 
(HZl) for C. 

Contribution III: We directly take the limit kq 
and we replace the factors B by their asymptotic expres- 
sions in (l?T|) . After the change of variable ka = /ciq, 
angular integration and averaging out of the oscillating 
terms fej*^^*" , this gives 



(niii(ki)) = 



9 4^5|_^^|2 



27r2 



kl 



+ 00 



(A9) 

In this result, we change the integration variable setting 
g = e", where a varies from —00 to +00. The odd com- 
ponent of the integrand (involving sinh a) gives a vanish- 
ing contribution. The even component of the integrand 
involves a rational fraction of cosh a to which we apply 
a partial fraction decomposition. Then we use (|A7p to 
obtain 



niii(fci 



47r5|A/'^|2 3V3sinh(7r|so|/3)-Fsinh(27r|so|/3) 



kl 



271 



sinh(7r|so|) 
+ o{l/kl). 



(AlO) 



Contribution IV: We directly take the limit ko — > 
and we replace the factors B by their asymptotic ex- 
pressions in (j22|) . We perform the change of variable 
k2 = (^i/2)q, we average out the oscillating terms kf^^° . 
The angular integration in spherical coordinates of axis 
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the direction of ki may be performed using 



6. Oscillating contribution at large fci 



/ 



dv 



1 



so/2 



(1 



1 



so/2 



/so, 



1 — V J \1 — V ^ 

(All) 

where the variable v is restricted to the interval (—1,1). 
This leads to 



(niv(ki)) = 



47r5|AAv,|2 36 



+ 00 



dq 



g + 1 



+ o{l/kl). (A12) 



Calculating this integral directly is not straightforward 
because of the occurrence of the absolute value \q — 1|. 
We thus split the integration domain in two intervals. 
For q E [0, 1] we set q — {X — 1)/{X + 1) (an increasing 
function of X, where X spans [1, +oo]). For q E [1, +oo] 
we set q = {X + 1)/{X — 1) (a decreasing function of X, 
where X here also spans [1, +oo]). Then 



(niv(ki)) 

(X^ - 1 + x-^ 
{X^ + l + X-^] 



k\ 2 



[s^'X^« 



^«]}+oil/kl). 

(A13) 



We then set X = e", where a ranges from zero to +oo, 
and we use the fact that the resulting integrand is an 
even function of a to extend the integral over the whole 
real axis. We integrate by parts, integrating the factor 
sin(a|so|), and we perform a partial fraction decompo- 
sition of the resulting rational fraction of cosh a. Using 
and its derivatives with respect to 6, we get 



(niv(ki)) = 



72^5 |_^^ 1 2 



sinh( 



27r|so| 



I cosh( 



27r|so| 



cosh( 



3 ' ' 3 ' 

)+sinh(^)]}/[2^|so|sinh(^|so|): 
+ o{l/kl). (A14) 



Sum of the four contributions: Summing up the non- 
oscillating terms in l/fcf of the contributions nj, nn, nm 
and niv, we obtain as a global prefactor 



5 = 



\/3 cosh(27r|so|/3) — cosh(7r|so|/3) 



(A15) 



|so| sinh(7r|so|) 

Multiplying ([3]) on both sides by sinh(|so|7r/2) and using 

2 sinha sinh & = cosh(a + 6) — cosh(a — 6), Va, 6 (A16) 

we find that S is exactly zero. As a consequence, the 
non-oscillating part of the momentum distribution of an 
infinite scattering length Efimov trimer behaves at large 
ki as 



(n(ki)) = ^ + o(l/fc^). 



(A17) 



In the large ki tail of the momentum distribution, we 
now include oscillating terms, having oscillating factors 
such as k^'^'^" . The calculation techniques are the same as 
in the previous subsection, so that we give here directly 
the result. We find that the leading oscillating terms 
scale as 1/fcf: 



i(ki)-(n(ki)) = -iy|A/-^P 



A 



kiVi 



2so 



-I- C.C. 



+ o{l/kl) (A18) 



where the complex amplitude A is the sum of the con- 
tributions coming from each of the four components 
P9I20I21I22P of the moment distribution. 



A^Ai 
We successively find 
3 



Au + Aui + Aiv. 



Ai = 



87r2 



+ 00 

1 



Wtt' 



3 /•+°°, q^'° 



Au = dq'i-^[iq' + q^ + ir 



(A19) 



(A20) 
1] (A21) 



V3 sinh(47r| So I /3) + sinh(27r I So I /3) 



A 



III 



47r 

3 f+°°, 
2^1 



sinh(27r|so| 



(A22) 



V3 



I sinh( 



q 

2^|so| 



) 4- sinh( 



cosh( 



27r|so| 



A 



IV 



+ 00 



dq 



1 47rsinh(7r|so|) 
3 ' 

-cosh(^)]}, (A23) 
q{l + q^Y" 



(g2+ 3)2(^2 + 1) 



2<?/(l+<?^) ti _ y2Nso/2 

dv 



0.0243657158 - 0.069868025H. 



(A24) 
(A25) 



We have calculated analytically all these integrals, ex- 
cept for (|A24[) where the angular integration gives rise to 
the integral over v and thus to a difficult hypergeometric 
function. Finally 



A ~ 0.1022397786 - 0.1218775240?. 



Appendix B: Direct calculation of E{A) 



(A26) 



To calculate the cut-off dependent energy E{A) defined 
in p2p for an infinite scattering length Efimov trimer, 
the method consisting in calculating the momentum dis- 
tribution n(k) and then integrating ([32]) is numerically 
demanding: A double integral has to be performed to ob- 
tain n(k), see (HH), so that the evaluation of E{A) results 
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in a triple integral. A more direct formulation, involving 
only a double integration, is proposed here. One simply 
rewrites (l32l) as 



E{A) 



(2^ 



f{k) 



2m 



n(k). 



C_ 

fc4 



(Bl) 



where the function /(fc) is equal to unity for < fc < A 
and is equal to zero otherwise. Then one plugs in (jBl[) the 
expression ([T5| of n(k), also replacing C with its integral 
expression (jA4l) . An integration over two vectors in M"^ 
appears, so that E{h) — £'oasy(A) + £^hard(A) with 



£'casy(A) 



3^ 

43 



(27r)3''^^ 2m 



(27r)3 



2B2(^g) 



(B2) 



•?4 /• j3 L 



2m 



_cfq_ 
(27r)3 



2^(^|q + k/2|)^(^|q-k/2|) 

(q2 + 3p ^ ^2)2 



(B3) 



The first part -Eeasy of this expression originates from the 
bits ni, nil, nm of the momentum distribution and from 
C; angular integrations may be performed, one is left 
with a double integral over the moduli k and q. Taking 
Kq as a unit of momentum and h^n^/m as a unit of energy 
in what follows: 



easy (A) — ( 



A U 



+ 00 



dqq 



2B\^q) 



fc4 



B^^k) + 252(^(7) + iBi^k)B{%) 



that we integrate numerically. The second part -Ehard (A) 
in (IB3p originates from the bit riiy of the momentum 
distribution. Performing the change of variables q = 
(ki — k2)/2 and k = ki + k2 ensures that the factors B 
are now functions of the moduli ki and ^2 only. 



-EhardCA) 



43 



d^ki f d^k2 1 



(2.)3y (2.)3 2(^^+'^^)^^(l'^^+^^l) 



2B{^ki)B{^k2) 



1 



(B5) 



[ki + /cl + ki • k2 

so that angular integrations may again be performed, in- 
volving the integral 

du ■ 

-1 V^i ^ "-2 



I^kiM) 



kl 



2kik2U 



{ki + A:| + kik2U - 



2kik2U 



1 



fclfc2 

H 



ln(l + fc^ + fc^ 



1) 



kik2u) 



(B6) 



fc|)/2 



i + ki 



ko + kik2U 



inax[— l,min(l,C/)] 



(B7) 



where u is the cosine of the angle between the vectors 
kl and kj, U = [A^ - {kf + fc|)]/(2fcifc2), max(a, &) 
(resp. min(a, 6)) is the largest (resp. smallest) of the 
two numbers a and b, and the notation stands 
for F{b) — F{a) for any function F{u). We also used the 
fact that |ki + k2| < A if and only ii u < U. This leads 
to 



-EhardCA) — ( — 



+00 



+00 



dki kl / dfc2 fc2 /(fci, fc2) 



xS(^fci)S(^fc2). (B8) 

Further simplifications may be performed. One can map 
the integration to the domain fci > fc2 since the integrand 
is a symmetric function of fci and fc2. Then performing 
the change of variable ki — q + fc/2 and k2 = q — k/2, 
and using the fact that /(fci, fc2) = if fci — fc2 > A, we 
obtain the useful form 



-Ehard(A) — ( — 

xI{q+h,q-h)B 



4 pA I'+oo 

2 dk dq (g2 - fc2/4)2 

Jo Jk/2 



^{<i + k/2) 



B 



%i~k/2) 



(B9) 



that we integrate numerically. A useful result to control 
the numerical error due to the truncation of the integral 
over g to a value ^ A and ^ 1 is I{q+k/2, q~k/2) ^ 

(fc4-A4)/(8g6). 



Appendix C: Analytical expression for E 

As explained in the main text, the naive regularization 
((33| of the energy formula gives an energy E that actually 
differs from the energy of the trimer -Etrinn because the 
momentum space cut-off function was introduced at the 
last stage of the calculation. Here we introduce a mo- 
mentum cut-off function at the level of the three-body 
wavefunction, simply by making the substitution 



B{k) B,,(k) = B{k), 



(CI) 



where we have set k — kq sinh a so that a = 
+ fc2/Kg -I- k/no]. For this consistent regulariza- 
tion, we expect that the usual energy formula holds in 
the limit of a vanishing ry, and this was checked explicitly 
in Appendix H of [l^. This means that 



lim E„ 

r)->-0+ 



Et, 



with 



E„ 



d^k h^e 



(27r)3 2m 



(C2) 



(C3) 
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The single-particle momentum distribution is ob- 
tained by replacing the function B by in 
Eas. (|18ll9l20l21l2"2|) . Its large k asymptotic behavior can 
be obtained along the lines of Appendix |XJ «,j(k) = 
rV^^y'^P{k) + 0{l/k^) with 



The last step is to calculate with the same tech- 
niques as in the Appendix El One finds 



27r2 



.4%/3 



(C7) 



asymp 



-2ri ln(y3fe/Ko) 



-\- D,j cos 



2 1 So I ln(\/3fc/Ko 



}. (C4) 



In the limit 77 — )■ 0+ one has to recover (|26p so that 
the coefficient Z)^ of the non-oscillating 1/fc^ contribution 
tends to zero in that limit. However, as we shall see, 
Dtj/t] does not tend to zero, which leads to the failure 
of the naive regularization. The expressions of the other 
coefficients C^, Z?,, and ipr/ are not needed. 

Following we define Snjj{k) as n^(k) — Cjj/k^ for 
k < /cmin and as n^(k) - n^''y'"P(fc) for fc > kmin- This 
results in the splitting 



(Pk h^e 



(27r)3 2m 

/.+00 



5n,i{k) 



47r2r 



da;e ^''^[A, + £',7 cos(2|so|x -I- .^^)] 

(C5) 



where the change of variable x = ln(-\/3fc/Ko) was used 
so that Xmin = ln('\/3fcmin/'«o)- For T] — >■ 0+, we can re- 
place in the right hand side of (IC5|) Snri(k) with 6n{k) 
since the first integral converges absolutely, but we can- 
not exchange the 77 — )■ 0+ limit and the integration in the 
second integral. After explicit calculation of this second 
integral, we take 77 0"*" and we recognize E from 
so that 



'trim 



E 



—5 — hm . 



with 



J 7^ 



-2?7 In q 



1 + 



+00 



The contri 



dq 



1 



l + q' 
8q 



-„l„(i^) 



l + g2 (g2 + 3)2 



_2q 

e 

aw — 



(C8) 
(C9) 



1 



1 



so/2 



C.C. 



(CIO) 



ibutions and iiT^ originate respectively 

from the bits njj , nm and n/\/ in the decomposition (|18p 
generalized to n^. Taking the derivative with respect to 
77 and then taking the limit 77 — > 0"*" gives More 
precisely, one finds that 



drj 
dJr, . 



+00 



dq{~\nq) 



1 + q^ +q 



(Cll) 
4 (C12) 



which, together with (|38p and (jC6|) , suffices to determine 
^^^|^_i.o+ so that we do not reproduce here its lengthy 

expression. The remarkable fact that --^\ri^o+ i^^^Y be 
expressed as an integral over a single variable q, whereas 
the expression of for a general rj in (jC10[) involves a 
double integral, results from an integration by part over q 



(C6) in (|C10p . taking the derivative of the bit 



29/(1+9') 



dv . 
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